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Abstract 

Global  kinetics  for  the  liquid  state  thermochemical  decomposition  of  neat  1,3, 5, 5- 
tetranitrohexahydropyrimidine  (DNNC),  perdeuterio-labeled  DNNC-d6,  and  partially  deuterium-labeled 
DNNC-d2  and  DNNC-d4  isotopomers  were  obtained  by  isothermal  differential  scanning  calorimetry  (IDSC). 
Molecular  kinetic  deuterium  isotope  effect  (KDIE)  values  obtained  with  DNNC  and  DNNC-d6  from  174  to 
194  °C  revealed  that  C-H  bond  rupture  regulates  both  an  endothermic  catalytic  initiation  and  the  exothermic 
propagation  of  the  liquid  thermochemical  decomposition  process.  Using  IDSC-based  KDIE  comparisons  with 
the  DNNC-d2,  DNNC-d4,  and  DNNC-d6  isotopomers,  a  more  detailed  chemical  structure/mechanistic 
relationship  emerged  by  entering  the  interior  of  the  DNNC  molecule.  Here  structural  kinetic  KDIE  results 
showed  the  rate-controlling  C-H  bond  rupture  has  its  origin  at  the  non-equivalent  C-2  methylene  group 
sandwiched  between  the  two  nitrated  DNNC  nitrogen  ring  atoms,  versus  at  the  chemically  equivalent  C-4  and 
C-6  methylene  ring  positions  located  elsewhere  in  the  DNNC  molecule.  Elucidation  of  such  mechanistic 
features  should  aid  in  the  structural  design  of  new  high  energy  compounds  with  improved  thermochemical 
properties.  A  170.0  kJ/mol  activation  energy  appeared  for  the  endothermic  induction  period,  and  a  lower 
104.2  kJ/mol  activation  energy  was  determined  for  the  exothermic  acceleratory  portion  of  the  DNNC 
decomposition  process.  The  global  liquid  and  solid  state  thermochemical  decomposition  processes  for  DNNC 
are  compared. 

Keywords:  1,3,5,5-Tetranitrohexahydropyrimidine;  DNNC;  Liquid  state  decomposition;  Deuterium  isotope 
effect;  KDIE 
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1.  Introduction 


Although  complex  global  reaction  kinetics,  obtained  by  relatively  inexpensive  thermal 
analysis  methods,  often  are  ignored  in  favor  of  highly-defined  individual  kinetic  steps 
requiring  more  sophisticated  and  expensive  equipment  [1],  global  kinetics  data  can  describe 
complex  thermochemical  decomposition  processes  in  ways  which  critically  impact  practical 
technology  and  safety  concerns  [1,  2].  For  example,  isothennal  differential  scanning 
calorimetry  (IDSC)  global  kinetics  rate  constants,  detennined  for  the  solid  state 
thermochemical  decomposition  of  1, 3, 5-triamino-2, 4, 6-trinitrobenzene  (TATB),  accurately 
predict  the  354  °C  experimentally  detennined  critical  temperature  which  initiates  a  thermal 
explosion  event  [3]. 


The  condensed  phase  neat  thermochemical  decomposition  process  an  energetic 
compound  follows  is  primarily  a  complex  mechanistic  set  of  interacting  intermolecular 
bimolecular  reactions  [4]  with  some  degree  of  intramolecular  unimolecular  character. 
Generally,  two  key  chemical  bond  ruptures  occur  during  the  thennochemical  decomposition 
process,  the  initial  bond  rupture  and  the  rate-controlling  bond  rupture;  very  rarely  does  one 
bond  rupture  constitute  both.  The  first  is  the  initial  homolytic  bond  rupture  (scission  or 
cleavage)  which _ 
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1  Portions  of  this  material  were  presented  at  the  following  conferences:  35th  International  Conference  of  ICT:  Energetic 
Materials,  Structure  and  Properties,  Karlsruhe,  Germany,  Jun  29-Jul  2,  2004;  203rd  Nat.  Am.  Chem.  Soc.  Mtg.,  San 
Francisco,  CA,  USA,  Apr  5-10,  1992;  and,  22nd  International  Conference  of  ICT:  Combustion  and  Reaction  Kinetics, 
Karlsruhe,  Germany,  Jul  2-5,  1991. 

2  Experimental  IDSC  thennochemical  decomposition  investigations  were  conducted  at  The  F.  J.  Seiler  Research 
Laboratory  (FJSRL),  USAF  Academy,  CO,  USA  prior  to  its  Sep  30,  1995  decommissioning;  (a)  Directorate  of  Chemical 
Sciences,  FJSRL,  (b)  Department  of  Chemistry,  USAF  Academy. 
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determines  the  specific  pathway  (e.g.  A,  B,  or  C  in  Fig.  1)  the  decomposition  process  will 
follow  and  the  resultant  product  composition.  Certain  structural  arrangements  in  an 
energetic  molecule  {e.g.  N-NO2  for  secondary  nitramines)  [5]  are  reaction  centers 
(explosophores)  [6]  for  this  initial  bond  rupture  which  often  are  termed  the  “trigger 
linkage,”  or  more  recently,  the  “first  reaction”  [7].  Once  this  specific  pathway  is  established 
(pathway  C  in  Fig.  1),  a  myriad  of  subsequent  chemical  reaction  channels  proceed  which 
occur  in  series  and  competing  parallel  sequences.  Along  this  reaction  pathway,  one 
individual  chemical  reaction  in  a  specific  reaction  channel  proceeds  more  slowly  than  all 
others  (longest  arrow)  where  individual  reaction  steps  are  represented  as  arrows  (Fig.  1).  It 
is  during  this  slowest  reaction  step  that  the  second  critical  homolytic  bond  rupture  occurs 
which  significantly  regulates  the  overall  rate  at  which  this  slowest  reaction  channel  and  the 
entire  decomposition  process  proceeds.  This  subsequent  KDIE-determined  rate-controlling 
chemical  bond  rupture,  accompanied  by  a  concomitant  atom  transfer  [8],  occurs  early  in 
the  thermochemical  decomposition  process  [2,9,10]  and  usually  is  not  synonymous  with  the 
trigger  linkage,  but  closely  follows  it.  This  rate-controlling  bond  rupture  can  be  determined 
in-situ  using  the  kinetic  deuterium  isotope  effect  (KDIE)  approach  in  conjunction  with 
ambient  pressure  isothermal  DSC  [2, 3, 8-1 1]  reaction  quenching  [1 1],  or  TGA  [12]  analyses. 
If  this  slowest  reaction  did  not  exert  a  significant  impact  on  the  overall  global 
decomposition  rate,  a  noticeable  global  KDIE  would  not  be  observed.  In  the  example 
illustrated  by  Fig.  1,  note  that  the  slow  individual  rate-controlling  reaction  step  rate  for 
either  the  normal  energetic  compound  or  its  deuterium  labeled  analogue  cannot  be  measured 
in  isolation  where  a  specific  KDIE  =  2.00  would  be  determined.  Rather,  the  entire  global 
decomposition  process  must  be  analyzed  and  the  overall  rates  for  the  energetic  compound 
and  its  deuterated  analogue  are  compared  when  determining  a  KDIE  value  [13].  Thus,  the 
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measured  global  KDIE  value  ( 1 .45)  is  likely  a  smaller  number  than  what  would  be  found  for 
the  specific  KDIE  value  (2.00)  in  the  individual  rate-controlling  step. 

Insert  Fig.  1 

This  critical  rate-controlling  bond  rupture  has  been  determined  for  the  ambient  pressure 
thermochemical  decomposition  processes  of  TNT  [10],  TATB  [3],  HMX  [8,12],  and  RDX 
[1 1,12].  In  every  case,  the  same  mechanistic  rate-controlling  bond  rupture  found  in  the 
ambient  pressure  decomposition  process  for  each  compound,  proved  to  be  a  significant  rate¬ 
controlling  mechanistic  feature  observed  in  rapid  deflagration  [8,14],  high-pressure 
combustion  [14-16],  and  thennal  explosion  events  [3,14].  Furthermore,  it  is  detected  as 
being  a  significant  contributing  mechanistic  feature  in  impact  or  shock  explosive  initiation 
sensitivity  12,17-23].  Identification  of  the  pathway  directing  initial  bond  rupture,  of  the 
subsequent  rate-controlling  bond  rupture,  and  of  the  resultant  product  composition  formed 
during  the  thermochemical  decomposition  processes  are  critical  for  obtaining  energetic 
compounds  with  improved  thermochemical  and  initiation  sensitivity  properties. 

The  first  KDIE  approach  applied  to  the  thermochemical  decomposition  of  an  energetic 
compound  occurred  with  neat  liquid  TNT  and  TNT-<2-d3  [10].  IDSC  analysis  during  the 
TNT  endothermic  induction  period,  the  earliest  portion  in  the  thermochemical 
decomposition  process,  revealed  that  C-H  bond  rupture  in  the  pendant  TNT  methyl  group 
controls  the  formation  of  a  catalytic  species,  which  upon  reaching  a  minute  threshold 
concentration,  initiates  the  TNT  exothennic  acceleratory  phase  and  a  concomitant  energy 
release.  A  subsequent  KDIE  study  with  solid  state  TATB  and  TATB-d6  found  that  N-H 
bond  rupture  also  constituted  the  rate-controlling  mechanistic  feature  early  during  the 
thermochemical  decomposition  process  [3]. 
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The  seminal  KD IE-based  thermochemical  decomposition  investigation  of  a  heterocyclic 
nitramine  was  accomplished  with  HMX  and  HMX-dx  and  revealed  a  different  KDIE  and 
rate-controlling  mechanistic  step  dependence  on  the  HMX  physical  state  [8].  This  IDSC 
study  showed  a  rate-controlling  C-H  bond  rupture  occurred  during  solid  state  decomposition 
along  with  other  rate-controlling  features  in  the  mixed-melt  and  liquid  state  decompositions. 
Proposed  was  an  initial  N-NO2  bond  rupture  to  fonn  the  NO2  radical  which  then,  either, 
unimolecularly  from  an  HMX  radical  intermediate  existing  in  a  proposed  cage  effect  [24],  or 
bimolecularly  from  a  neutral  HMX  molecule,  transfers  a  hydrogen  atom  via  NO2  radical 
abstraction  to  form  the  unstable  HONO  specie  [8].  This  unstable  species  undergoes  further 
decomposition  to  fonn  other  oxidizing  and  catalytic  species  [8],  such  as  NO,  NO2,  and  H20, 
or  in  some  cases,  possibly  OH  and  NO  [25],  while  resultant  HMX  radicals  continue  along 
various  decomposition  reaction  channels  and  condensed  phases  [8], 

This  rate-controlling  C-H  bond  rupture  for  solid  state  HMX  was  confirmed  by  a 
subsequent  HMX/HMX-dx  isothennal  TGA  decomposition  study  [12].  Later,  a  more 
detailed  mechanistic  isothennal  pyrolysis  investigation  with  HMX  and  HMX-dg  was 
conducted  using  simultaneous  thennogravimetric  modulated  beam  mass  spectrometry 
(STMBMS)  where  different  reaction  channels  of  specific  product  specie  formation,  during 
the  decomposition  process,  could  be  measured  concunently  [26,27].  At  least  four  different 
reaction  channels  were  identified,  and  a  large  KDIE  was  found  for  water  product  formation. 
Finding  a  large  KDIE  for  water  formation  also  was  consistent  with  an  initial  N-NO2  bond 
cleavage,  a  subsequent  rate-controlling  methylene  C-H  bond  rupture  and  hydrogen  atom 
transfer  from  an  HMX  species  to  form  the  unstable  HONO  molecule  [8]  and  its 
decomposition  to  the  water  byproduct.  Water,  occluded  in  the  HMX  crystal  structure,  also 
was  suggested  as  a  hydrogen  atom  donating  catalyst  for  another  reaction  channel  [27]. 


Distribution  A.  Approved  for  public  release;  distribution  unlimited. 


5 


Conceivably,  the  water  product,  fonned  via  the  KDIE  measured  rate-controlling  C-H  bond 
rupture  outlined  above,  also  could  contribute  to  the  water  catalysis  channel. 

The  same  rate-controlling  C-H  bond  rupture  was  identified  for  the  liquid  state 
thermochemical  decomposition  of  RDX  and  RDX-d6  using  isothermal  TGA  [12]  This  result 
was  further  verified  by  a  subsequent  IDSC  and  a  UV-based  reaction  quenching  KDIE  study 
[11].  A  later  STMBMS  RDX/RDX-d6  investigation  revealed  a  complex  multi-channeled 
liquid  thermochemical  decomposition  process  which  was  different  from  HMX,  but  still 
revealed  a  rate-controlling  C-H  bond  rupture,  and  atom  transfer  from  the  RDX  methylene 
group,  both  being  highly  prevalent  mechanistic  features  [25,28,29].  In  fact,  it  appears  the 
progression  rate  of  all  four  reaction  channels  may  depend  upon  C-H  bond  rupture  in  the 
first  channel. 

While  much  more  thermally  stable  than  the  RDX,  the  cyclic  DNNC  nitramine  derivative 
(Fig.  2)  undergoes  extensive  solid  state  thermochemical  decomposition  during  long-tenn 
exposure  in  the  142  to  148  °C  temperature  range  [2].  IDSC-based  KDIE  results  suggest  that 
ring  C-N  bond  rupture  for  the  solid  state  decomposition  likely  is  the  rate-controlling  bond 
rupture  during  the  long  endothermic  induction  period,  wherein  a  catalytic  species  likely 
forms,  and  upon  reaching  a  threshold  concentration,  initiates  an  autocatalytic  exothermic 
acceleratory  phase.  Once  the  acceleratory  phase  is  initiated,  C-H  bond  rupture  is  the  rate¬ 
controlling  mechanistic  feature  which  regulates  its  energy  release  rate.  The  endothermic 
induction  period  for  the  solid  state  decomposition  process  has  a  higher  activation  energy  of 
128  kJ/mole  compared  to  the  115-106  kJ/mole  [2,30]  determined  for  the  exothermic 
acceleratory  phase.  This  supports  the  likely  formation  of  a  catalytic  species  from  the  DNNC 
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molecule  during  the  solid  state  induction  period  and  mirrors  the  thermochemical  behavior 
followed  by  the  neat  liquid  TNT  thermochemical  decomposition  process  [10]. 


The  relationship  of  the  mechanistic  rate-controlling  homolytic  bond  rupture  found  in 
ambient  pressure  thermochemical  decomposition  processes  to  higher  temperature,  higher 
pressure  catastrophic  events  such  as  deflagration,  combustion,  thermal  explosion, 
detonation,  and  to  energetic  compound  sensitivity  issues,  is  well  documented  [9,14].  Given 
the  importance  of  this  high  energy  rate-controlling  mechanistic  feature,  the  thermochemical 
decomposition  of  neat  liquid  DNNC,  and  its  three  isotopomers,  DNNC-d6,  DNNC-d/t,  and 
DNNC-di,  was  investigated  by  isothennal  DSC  to  determine  the  rate-controlling  homolytic 
bond  rupture  and  its  location  in  the  DNNC  molecular  structure.  Similarities  and  differences 
between  this  neat  liquid  state  DNNC  thermochemical  decomposition  process,  and  the 
previously  reported  mechanistic  solid  state  DNNC  thermochemical  decomposition 
investigation  [2],  also  are  compared. 


2.  Experimental 

Caution!  Although  reasonably  stable  to  accidental  explosive  initiation,  DNNC  is  a  high 
energy  explosive  material.  Care  should  be  taken  to  use  proper  laboratory  shielding  and 
safety  procedures  in  all  handling  and  transfer  operations. 

2. 1  Isothermal  DSC  analysis 

All  isothermal  DSC  measurements  were  made  with  a  Perkin-Ehner  DSC  7  instrument 
described  previously4  using  DNNC  sample  sizes  ranging  from  1.50  to  1.53  mg.  Prior  to 
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conducting  the  thennochemical  decomposition  runs,  the  DSC  7  temperature  was  calibrated 


at  a  1  deg  K/minute  scanning  rate  using  a  two-temperature  calibration  procedure  with 
indium  (mp.  156.60  °C)  and  tin  (231.88  °C)  standards.  The  temperature  of  the  DNNC 
sample  (mp  =  154.8-156.2  °C)  [31]  was  stabilized  to  30  °C  under  the  nitrogen  gas 
atmosphere  and  then  was  brought  to  its  pre-determined  isothermal  decomposition 
temperature  (174,  178,  182,  186,  190,  or  194  °C)  at  a  480  °C  per  minute  heating  rate. 

2.2  Isothermal  DSC  data  reduction 

The  length  of  the  endothennic  induction  period,  and  resultant  induction  times,  were 
visually  determined  from  individual  thermograms  generated  for  each  IDSC  run.  The  minute 
units  were  then  converted  to  seconds  for  the  calculations  used  to  obtain  activation  energy 
and  KDIE  values.  For  determining  the  exothermic  acceleratory  phase  rate  constants, 
software  programs  were  written  at  our  laboratory  to  transfer  the  data  from  each  designated 
IDSC  run  to  a  separate  PC  station,  and  then,  to  input  the  data  into  an  autocatalytic  kinetics 
plot  program  that  used  the  5a/5 1  =  ka(l-a)  autocatalytic  rate  equation,  where  a  is  the  mole 
fraction  of  decomposed  DNNC  at  any  given  time  (t).  A  software  program,  based  on  that 
developed  at  the  Los  Alamos  National  Laboratory  [32],  was  used  and  is  described  in  detail 
elsewhere  [2,8]. 

2.3  Synthesis  of  DNNC  and  the  DNNC-dr„  DNNC-d4,  and  DNNC-dj  structural  isotopomers 
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The  synthesis  of  ADIOL  and  ADIOL-d4  (2,2-dinitropropane-l,3-diol-l,l,3,3-d4) 
reactants,  needed  to  obtain  DNNC  and  its  DNNC-d6  analogue,  was  previously  described 
[33].  The  partially  deuterium-labeled  DNNC-d2  and  DNNC-d4  compounds  were 
synthesized  using  an  identical  Mannich  condensation  reaction  procedure,  and  scaled  as 
needed,  by  employing  the  appropriate  choice  of  unlabeled  or  deuterium  labeled  2,2- 
dinitropropane-l,3-diol  (ADIOL  and  ADIOL-d4)  and  fonnaldehyde  or  formaldehyde^, 
plus  unlabeled  /er/-butylamine  in  methanol  solvent.  These  initial  Mannich  condensation 
reactions  provided  the  DBNP-d2  and  DBNP-d4  intermediate  compounds  which  were  then 
converted  to  their  respective  final  DNNC-d2  and  DNNC-d4  products  by  100%  anhydrous 
nitric  acid  nitration.  All  ’H  NMR  (acetone-d6,  TMS  std.),  mass  spectral,  FTIR 
(Kubelka/Munk)  powder  analyses  and  melting  points  for  DBNP-d2,  DBNP-d4,  DNNC-d2 
(155-157  °C),  and  DNNC-d4(155-156  °C)  were  consistent  when  compared  with  previously 
reported  values  for  the  unlabeled  DBNP,  and  DNNC  analogues  [33].  Yields  obtained  for 
DBNP-d2  and  DBNP-d4  compounds  were  similar  to  DBNP  (75%)  and  DBNP-d6  (72%). 

Final  product  yields  for  DNNC-d2  (63%)  and  DNNC-d4  (75%)  compared  well  with  DNNC 
(72%)  and  DNNC-d6  (71%)  [33].  Within  the  accuracy  and  precision  of  !H  NMR  integration 
values,  the  intermediate  DBNP-d2  reflected  an  incorporation  of  2.2  methylene  deuterium 
atoms;  DBNP-d4  incorporated  4.0  methylene  deuterium  atoms;  and,  5.9  methylene 
deuterium  atoms  were  incorporated  into  the  DBNP-d6  molecule  relative  to  the  two  unlabeled 
tert- butyl  groups  substituted  on  the  ring  nitrogen  atoms  at  the  1-  and  3-ring  positions. 
Because  of  the  synthesis  method  used  [33],  the  -CD2-  labeled  methylenes  were  incorporated 
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into  the  DNNC  ring  structure  as  a  discreet  unit;  mixed  -CHD-  methylenes  could  not  fonn. 


Proton  NMR  integrations  of  the  two  different  methylene  groups  (5  6.30  C-2  and  6  5.47  C-4 
&  C-6)  for  both  DNNC-d4  and  DNNC-d2  remained  in  a  1:2  integration  ratio,  respectively, 
indicating  a  random  distribution  of  the  labeled  -CD2-  deuterium-labeled  methylene  groups 
occurred  at  all  three  C-2,  C-4,  and  C-6  ring  positions.  In  spite  of  this  randomized  -CD2- 
methylene  distribution,  relevant  mechanistic  information  was  obtained  by  comparing 
DNNC,  DNNC-d2,  DNNC-d4,  and  DNNC-d6  IDSC  isotopomer  induction  times  and  rate 
constants. 

3,  Results  and  discussion 

The  initial  homolytic  bond  rupture  of  an  energetic  material  controls  the  reaction 
pathway,  the  number  of  reaction  channels,  the  decomposition  products,  and  the  product 
composition  of  a  thermochemical  decomposition  process.  Later  rate-controlling  homolytic 
bond  rupture  ultimately  controls  the  overall  energy  release  rate  during  the  decomposition 
of  a  given  energetic  material.  Improved  thermochemical  properties  must  start  with 
decreasing  the  rate  at  which  a  decomposition  process  proceeds,  either  kinetically,  or  by 
chemical  structure  alteration.  An  excellent  example  of  possible  progress  is  found  in  the 
chemical  structure  similarity  and  differences  which  exist  between  the  symmetrical  RDX 
cyclic  nitramine  molecule  and  the  unsymmetrical  DNNC  cyclic  pseudo-nitramine  (Fig.  2). 
This  led  us  to  investigate  the  solid  and  liquid  state  decomposition  of  1, 3,5,5- 
tetranitrohexahydropyrimidine  (DNNC)  pseudo-nitramine  [2,30]  for  comparison  with  RDX. 

First  synthesized  in  1984  [34-36],  the  DNNC  molecule,  is  a  unique  heterocyclic  hybrid 
structure  which  contains  structural  features  found  in  both  saturated  cyclic  nitramines  (RDX 
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and  HMX)  and  in  aliphatic  geminal-dinitroalkyl  compounds  (Fig.  2).  While  RDX  contains 
three  chemically  equivalent  (-CH2-)  methylenes,  the  DNNC  structural  variation  provides 
three  methylene  (-CH2-)  groups  where  only  two  are  chemically  equivalent.  Of  the  three 
DNNC  methylenes  situated  at  the  C-2,  C-4,  and  C-6  positions,  the  two  methylene  groups  at 
the  C-4,  and  C-6  positions,  adjacent  to  the  geminal  dinitro-substituted  C-5  carbon  [- 
C(N02)2-],  are  in  a  chemically  equivalent  environment  like  that  found  in  RDX,  but  the  third 
DNNC  methylene  group  at  the  C-2  position,  sandwiched  between  the  two  nitrated  ring  N 
atoms  (=NN02),  is  chemically  non-equivalent  with  respect  to  the  C-4  and  C-6  methylene 
positions. 

Insert  Fig.  2 

Because  of  this  structural  variation,  DNNC  thermochemical  stability  resembles  that  of 
TNT,  rather  than  its  structurally  related  RDX  and  HMX  nitramines,  but  does  provide  an 
explosive  energy  output  which  rivals  RDX  [2,33].  Like  TNT  (mp  =  80.0  to  81.2  °C)  [10], 
neat  DNNC  melts  without  noticeable  decomposition  when  heated  above  its  154.8  to  156.2 
°C  melting  point;  thus,  a  5.19  kcal/mole  (21.72  kJ/mol)  DNNC  heat  of  fusion  could  be 
detennined  by  DSC  analysis  [31]  (TNT  AHfusion  =  5.34  kcaFmol  or  22.34  kJ/mol)  [37]. 
Related  cyclic  nitramines  RDX  and  HMX  melt,  or  undergo  liquefaction,  with  decomposition 
rapidly  enough  that  their  heat  of  fusion  cannot  be  obtained.  Additionally,  DNNC  and  RDX 
structural  variation  affects  thermochemical  properties  related  to  rapid  pyrolysis  and 
simulated  combustion  behavior  [4,38]. 

The  liquid  state  thennochemical  decomposition  process  for  both  DNNC  and  DNNC-d6 
was  conducted  at  4  °C  temperature  intervals  from  174  °C  to  194  °C  using  isothermal 
differential  scanning  calorimetry  (IDSC).  Each  DNNC  and  DNNC-d6  IDSC  run  provided  a 
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thermogram  curve  displaying  three  main  features:  (1)  an  endothermic  induction  period,  (2) 
followed  by  the  initiation  of  an  exothermic  acceleratory  phase  which,  after  reaching  a 
maximum  energy  release  at  the  thermogram  peak,  (3)  entered  into  a  decreasing  exothermic 
decay  phase  until  complete  decomposition  resulted.  Ample  thermogram  curve  shifts 
between  each  IDSC  thennogram  were  obtained  at  the  six  isothennal  temperatures  (174,  178, 
182,  186,  190,  and  194  °C)  selected  for  relevant  kinetic  evaluation  (Fig.  3). 

Insert  Fig.  3 

3.1  Thermochemical  physical  and  kinetic  properties 

Within  this  20  °C  temperature  regime  (174  -  194  °C),  post  IDSC  sample  weights  of 
decomposed  neat  liquid  DNNC  showed  an  average  81%  mass  loss  from  conversion  to 
gaseous  products  (Table  1).  The  remaining  19%  of  the  DNNC  sample  was  converted  to 
solid  phase  residue  product(s)  and  remained  in  the  sample  pan.  The  liquid  DNNC-d6 
decomposition  process  gave  a  nearly  identical  weight  loss  average  of  80%  with  a  20%  solid 
phase  residue  (Table  2).  For  comparison,  solid  state  DNNC  thennochemical  decomposition 
afforded  a  76%  mass  loss  and  solid  state  DNNC-d6  a  slightly  smaller  73%  mass  decrease. 
Formation  of  a  solid  phase  residue  has  very  important  implications  in  IDSC  kinetic  data 
evaluation. 

Post  solid  phase  residue  can  mean  that  intermediate  derivatives  form  during  the 
thennochemical  decomposition  process  which  react  further  until  decomposition  is  complete. 
Two  problems  can  result,  and  both  can  cause  a  misinterpretation  of  the  KDIE  value  and  its 
associated  kinetics.  First,  kinetic  evaluation  may  provide  a  combined  rate  constant  value  to 
which  decomposition  of  the  subject  compound  and  further  decomposition  of  the 
intennediate  condensed  phase  derivatives  both  contribute.  This  results  in  a  KDIE  value 
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which  is  diluted  or  reduced  from  the  true  value  of  the  subject  compound.  Secondly,  this  rate 
constant  may  reflect  data  from  decomposition  of  the  condensed  phase  derivatives  only 
formed  from  the  initial  energetic  compound.  In  this  case,  the  KDIE  value  offers  no  useful 
information  about  the  original  energetic  compound  and  its  decomposition  process.  These 
two  related  problems  are  amply  documented  in  ID  SC  thermochemical  decomposition 
investigations  with  liquid  TNT  [9,10]  and  solid  TATB  [3,9].  When  the  presence  of 
condensed  phase  products  appear,  only  the  early  portion  of  the  thermochemical 
decomposition  process,  and  the  corresponding  thermogram  portion,  should  be  evaluated, 
notably,  the  endothermic  induction  period  and  the  exothennic  acceleratory  phase. 

The  DNNC  and  DNNC-d6  thermochemical  decomposition  process  afforded  smooth 
thennograms  that  followed  an  autocatalytic  rate  behavior  during  a  large  portion  of  the 
exothermic  acceleratory  phase  (Fig.  3).  The  obvious  incremental  rate  increase  as  one 
proceeds  to  higher  temperatures,  is  visually  depicted  with  decreasing  induction  times  (t,)  and 
steeper  (faster)  acceleratory  phase  portions  reflecting  a  higher  rate  constant  (k).  The  percent 
of  the  acceleratory  phase  that  could  be  fitted  linearly  to  the  autocatalytic  rate  plot,  the 
exceptional  least  squares  linear  curve  fits  (R  )  obtained,  and  the  corresponding  rate 
constants  (k)  determined  from  the  slopes  of  the  linear  portion  of  the  autocatalytic  rate  plot  at 
each  temperature,  are  found  in  Table  1.  The  rate  constants  shown  at  each  isothermal 
temperature  are  averaged  over  at  least  five  runs.  Fig.  4  illustrates  the  conversion  of  one 
thermogram  run  at  182  °C  into  its  autocatalytic  rate  plot.  Analogous  data  for  DNNC-d6 
appear  in  Table  2.  At  a  liberal  99%  Student  t  variation  analysis,  each  rate  constant  is 
unique  and  does  not  overlap  with  ones  at  adjacent  temperatures.  The  heat  of  decomposition 
for  both  DNNC  and  DNNC-d6  give  similar  values. 
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Insert  Fig.  4 


Compared  to  the  solid  state  DNNC  thermochemical  decomposition  process  [2],  both  the 
induction  period  and  acceleratory  phase  kinetics  proceed  much  faster  above  the  154.8  to 
156.2  °C  DNNC  melting  point.  For  example,  the  184  min  (1 1,040  s)  induction  time  for 
solid  state  DNNC  at  148  °C  is  nearly  7.5  times  longer  than  the  148  s  for  liquid  state  DNNC 
at  182  °C.  For  the  liquid  state  DNNC  acceleratory  phase  at  182  °C,  the  rate  constant  seen  in 
Table  1  is  slightly  over  one  order  magnitude  faster  than  the  3.38  x  10"4  s'1  rate  constant 
obtained  at  148  °C  for  solid  state  DNNC  [2],  The  liquid  state  DNNC  thermochemical 
decomposition  releases  heat  energy  much  more  rapidly  than  does  the  solid  state  process. 


Table  1." 

Physical  and  kinetic  properties  of  neat  liquid  DNNC  thermochemical  decomposition 


Temp. 

(°C) 

Induction  time 

(s) 

Rate  constant 
(kx  10‘3  s'1) 

k  linearity 
(%)b  ' 

AH  decomp. 

(kJ/mol) 

Wt.  loss 
(%) 

No. 

runs 

194 

52.6  ±28.3 

12.15  ±0.48 

68 

-802.5  ±48.5 

80  ±2 

8 

190 

72.2  ±  19.3 

9.61  ±0.43 

83 

-775.3  ±56.9 

82  ±4 

5 

186 

93.8  ±  17.5 

7.55  ±0.64 

94 

-801.6  ±  30.1 

82  ±4 

5 

182 

148.4  ±  13.8 

5.87  ±0.39 

95 

-786.6  ±38.1 

80  ±5 

5 

178 

228.6  ±  76.4 

4.58  ±0.25 

94 

-799.1  ±61.1 

80  ±7 

5 

174 

359.6  ±  153.7 

3.45  ±0.16 

95 

-803.3  ±38.1 

80  ±3 

5 

ave. 

-795.0  ±45.6 

81  ±4 

a  All  data  variances  computed  at  the  99%  Student  t  level. 

b  Percent  of  thermogram  exothermic  acceleratory  phase  fitting  linear  portion  of  the 
autocatalytic  rate  plot  at  a  0.999  least  squares  value  (R2) 
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Table  2a 

Physical  and  kinetic  properties  of  neat  liquid  DNNC-d6  thermochemical  decomposition 


Temp. 

(°C) 

Induction  time 

(s) 

Rate 

constant(k  x 
10-3  s'1) 

k  linearity 
(%)*  ' 

AH  decomp 
(kJ/mol) 

Wt.  loss 
(%) 

No. 

runs 

194 

73.6  ±35.1 

7.16  ±0.89 

57 

-766.5  ±  36.4 

80  ±3 

6 

190 

126.5  ±28.3 

5.94  ±0.45 

72 

-797.9  ±  119.7 

84  ±  12 

5 

186 

169.9  ±39.5 

4.71  ±0.35 

81 

-777.8  ±81.2 

81  ±4 

6 

182 

274.0  ±59.6 

3.58  ±  0.19 

82 

-767.8  ±53.3 

81  ±6 

5 

178 

437.4  ±  120.6 

2.76  ±0.29 

82 

-754.8  ±25.1 

81  ±4 

5 

172 

652.0  ±82.2 

2.06  ±0.10 

83 

-769.0  ±  24.3 

77  ±  1 

5 

ave. 

value 

-772.4  ±  56.9 

80  ±5 

— 

a  All  data  variances  computed  at  the  99  %  Student  t  level. 

b  Percent  of  thermogram  exothermic  acceleratory  phase  fitting  linear  portion  of  the 
autocatalytic  rate  plot  at  a  0.999  least  squares  value  (R2). 


Activation  energies,  computed  at  a  95%  Excel1  standard  error  level,  reveal  a  170.0  ± 
27.7  kJ/mol  energy  of  activation  (R  =  0.9922)  for  the  DNNC  induction  period  using  the 
five  lower  temperatures  from  174  to  190  °C  (Fig.  5).  At  higher  temperatures,  the  induction 
period  time  is  short  enough  that  it  approaches  the  instrumental  temperature  equilibration 
regime,  and  accurate  measurement  becomes  difficult.  Little  time  is  left  for  a  linear 
induction  period  portion  in  the  thermogram  curve  to  develop,  and  this  can  mask  the  point  of 
curve  deflection  where  the  acceleratory  phase  is  initiated.  An  induction  period  activation 
energy,  which  includes  all  six  temperatures,  equals  160.5  kJ/mol  but  gives  a  less  desirable 
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2 

linear  curve  fit  (R  =  0.9903);  therefore,  the  higher  170.0  kJ/mol  is  selected  as  more  valid. 
The  170  kJ/mol  graphic  plot  uses,  as  the  shortest  induction  time,  a  72.2  s  value  determined 
at  190  °C  which  is  in  the  same  time  regime  as  that  of  the  shortest  DNNC-d6  induction  time 
value  (73.6  s)  at  the  higher  194  °C  temperature. 

Insert  Fig.  5 

A  noticeably  lower  104.2  ±3.9  kJ/mol  energy  of  activation  is  found  for  the  DNNC 
exothermic  acceleratory  phase.  Using  the  rate  constant  values  for  all  six  temperatures,  an 
excellent  linear  data  fit  results  (R  =  0.9993)  as  seen  in  Fig.  6.  The  significantly  lower  65.8 
kJ/mol  energy  of  activation  (170.0  kJ/mol  -  104.2  kJ/mol)  for  the  exothermic  acceleratory 
phase  follows  the  same  trend  as  that  seen  for  the  solid  state  DNNC  thermochemical 
decomposition  process  [2].  This  trend  suggests  the  possible  formation  of  a  small  amount  of 
decomposition  product  from  DNNC  itself  during  the  endothermic  induction  period,  which, 
upon  reaching  a  threshold  concentration,  catalyzes  initiation  of  the  autocatalytic  exothennic 
acceleratory  phase  where  DNNC  releases  its  stored  energy.  This  lower  energy  of  activation 
and  autocatalytic  decomposition  process  for  the  exothermic  acceleratory  phase  of  liquid 
state  DNNC  fits  with  catalytic  behavior.  It  also  mirrors  the  same  pattern  as  that  seen  for  the 
solid  state  DNNC  decomposition  process  [2],  as  well  as  that  reported  for  the  neat  liquid 
TNT  catalyzed  decomposition  process  [10], 

Insert  Fig.  6 

Because  significantly  slower  kinetics  result  with  DNNC-d6,  induction  times  are  longer  at 
corresponding  temperatures  used  with  DNNC;  so,  all  six  temperatures  are  included  to  obtain 
the  179.2  ±  14.3  kJ/mol  (R2  =  0.9967)  induction  period  activation  energy  (Fig.  5).  The  same 
activation  energy  trend  was  detennined  for  the  neat  liquid  state  DNNC-d6  thermochemical 
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decomposition  process  (Table  3).  The  DNNC-d6  rate  constants  at  corresponding 
temperatures  are  significantly  lower  than  those  for  DNNC  and  reflect  the  presence  of  a 
positive  kinetic  deuterium  isotope  effect  (KDIE).  As  with  the  DNNC  thermochemical 
decomposition  process,  a  large  portion  of  the  DNNC-d6  exothermic  acceleratory  phase 
provides  an  excellent  R  =  0.999  linear  fit  with  the  autocatalytic  rate  plot.  A  noticeably 
lower  104.5  ±  10.4  kJ/mol  activation  energy  (R~  =  0.9949)  also  was  observed  for  the 
DNNC-d6  acceleratory  phase  compared  to  its  induction  period  (Figure  6).  For  both  DNNC 
and  DNNC-d6,  the  higher  induction  period  activation  energy  is  statistically  significant  even 
at  a  most  strenuous  99%  variance  level  (Table  3). 

Activation  energy  (Eact)  comparisons  for  solid  and  liquid  state  DNNC  thermochemical 
decomposition  also  are  shown  in  Table  3.  While  DNNC  is  considered  to  be  a  thermally 
stable  compound,  given  enough  time  and  energy,  solid  state  DNNC  proceeds  through  a  slow 
thermochemical  decomposition  process;  and,  like  its  liquid  phase  decomposition  process, 
proceeds  with  an  induction  period  activation  energy  that  is  higher  than  that  seen  for  the 
acceleratory  phase. 


Table  3.fl 

Activation  energy  (kJ/mol)  comparisons  between  liquid  and  solid  state  DNNC 
thermochemical  decomposition  processes 

Thermogram  Liq.  DNNC  (R“)  Liq.  DNNC-d6  (R“)  Sol.  DNNC  (R2)  Sol.  DNNC-d6 

portion  (R2) 

Induction  period  170.0  ±  50.8  (0.992)  179.2  ±  23.7  (0.997)  128.1  (0.992)  [2]  not  determined 

Acceleratory  phase  104.2  ±  6.5  (0.999)  104.5  ±  17.2  (0.995)  105.6  (calc.)  [2]  1 15.2  (0.998)  [2] 

a  Variances  computed  at  Excel  99%  standard  error  level 

3.2  Molecular  Kinetic  Deuterium  Isotope  Effect 
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The  molecular  KDIE  for  a  given  compound  is  determined  by  comparing  thermochemical 
decomposition  rates  of  the  un-labeled  molecule  (DNNC)  with  the  perdeuterio-labeled 
molecule  (DNNC-dr,).  Comparing  both  induction  period  and  acceleratory  phase  data  at  the 
same  temperatures  in  Tables  1  and  2,  shows  that  the  non-isotopically  labeled  (normal) 
DNNC  thermochemically  decomposes  more  rapidly  than  its  fully  deuterium-labeled  DNNC- 
d6  analogue.  Like  that  previously  reported  for  liquid  TNT  [10],  solid  TATB  [3],  solid  and 
liquid  HMX  [8,12],  solid  and  liquid  RDX  [1 1,12],  and  solid  DNNC  [2],  a  positive  molecular 
kinetic  deuterium  isotope  effect  (KDIE)  defines  the  liquid  state  DNNC  thermochemical 
decomposition  process. 

Insert  Fig.  7 

This  positive  KDIE  clearly  is  observed  in  Fig.  7  where  individual  DNNC  and  DNNC-d6 
thennograms  taken  at  174  °C  are  compared.  Fig.  8  displays  how  the  magnitude  of  the  KDIE 
can  reveal  the  specific  rate-controlling  homolytic  bond  cleavage  associated  with  the  DNNC 
molecule.  Table  4  lists  the  KDIE  values  obtained  at  each  temperature  using  the  DNNC  and 
DNNC-d6  induction  periods  and  the  overall  KDIE  determined  for  the  entire  172  to  194  °C 
temperature  range.  The  induction  period  KDIE  is  determined  by  dividing  the  DNNC-d6 
induction  time  with  the  DNNC  induction  time  (Wkh)  [9,10].  A  molecular  KDIE  value  of 
1 .76  found  during  the  induction  period  reflects  a  primary  effect  and  identifies  C-H  bond 
rupture  in  the  ring  methylene  (-CH2-)  groups  of  the  DNNC  molecule,  or  a  structural 
fragment  derived  from  it,  as  being  the  rate-controlling  mechanistic  feature  (Fig.  8). 

Insert  Fig.  8 

Table  4.a 

Induction  period  kinetic  deuterium  isotope  effect  (KDIE)  for  DNNC 
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Temperature 

(°C) 

DNNC-d6  induction  period 

(min) 

DNNC  induction  period 

(min) 

KDIE  value 

(no  units) 

194 

1.23  ±0.53 

0.88  ±0.17 

1.40  ±0.47 

190 

2.11  ±0.41 

1.20  ±0.29 

1.76  ±0.32 

186 

2.83  ±0.59 

1.56  ±0.41 

1.81  ±0.34 

182 

4.57  ±0.89 

2.47  ±0.10 

1.85  ±0.20 

178 

7.29  ±  1.79 

3.81  ±  1.13 

1.91  ±0.38 

172 

10.87  ±  1.21 

5.99  ±2.29 

1.85  ±0.30 

average 

1.76  ±0.33 

a  Data  variances  computed  at  the  99%  Student  t  level 


Immediately  following  the  endothennic  induction  period,  the  DNNC  thermochemical 
decomposition  process  enters  into  an  exothermic  autocatalytic  acceleratory  phase.  In  this 
portion  of  the  thermogram  (Table  5),  KDIE  behavior  is  determined  by  dividing  the  normal 
DNNC  rate  constant  by  the  deuterium-labeled  DNNC-d6  rate  constant  (kh/kd)  [10].  Within 
the  174  to  194  °C  temperature  range,  a  substantial  overall  molecular  KDIE  equal  to  1.65  is 
found  (Table  5).  Among  the  myriad  of  chemical  reaction  channels  in  play  during  this 
complex  exothennic  acceleratory  portion  of  the  thermochemical  decomposition  process,  a 
rate-controlling  C-H  bond  rupture  occurs  in  the  substantially  slowest  one.  This  C-H  bond 
rupture  is  derived  from  the  methylene  sites  of  the  cyclic  DNNC  molecule,  or  a  structural 
fragment  derived  from  it,  and  significantly  controls  the  overall  exothermic  energy  release 
rate.  This  1.65  KDIE  value  is  calculated  to  be  3.17  when  mathematically  normalized  to 
room  temperature  conditions  (298  K)  [39]  and  further  verifies  assignment  of  a  rate¬ 
controlling  C-H  bond  rupture  where,  at  298  K,  a  primary  KDIE  equal  to  2.5  is  considered  to 
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be  the  minimum  [40].  This  primary  1.65  KDIE  value  for  the  exothermic  acceleratory  phase 
decomposition  of  liquid  DNNC  is  quite  similar  to  that  found  for  the  liquid  decomposition 
process  of  structurally- 


Table  5." 

Acceleratory  phase  kinetic  deuterium  isotope  effect  (KDIE)  for  DNNC 


Temperature 

(°C) 

DNNC  rate  constant 

(kx  10'3  s'1) 

DNNC-d6  rate  constant 

(kx  10'3  s'1) 

KDIE 

expl.  value 

KDIE 

25  °C  [i 

194 

12.15  ±0.48 

7.16  ±0.89 

1.70  ±0.13 

3.36 

190 

9.61  ±0.43 

5.94  ±0.45 

1.62  ±0.09 

3.16 

186 

7.55  ±0.64 

4,71  ±0.35 

1.60  ±  0.1 1 

3.10 

182 

5.87  ±0.39 

3.58  ±0.19 

1.64  ±0.08 

3.14 

178 

4.58  ±0.25 

2.76  ±0.29 

1.66  ±0.12 

3.14 

174 

3.45  ±0.16 

2.06  ±0.10 

1.67  ±0.07 

3.12 

average 

1.65  ±0.09 

3.17 

a  All  data  variances  computed  at  the  99%  Student  t  level 


related  RDX  which  was  measured  using  three  different  analysis  methods  (Table  6). 
Rounded  off  to  the  nearest  tenth  decimal  place,  a  1 .6  molecular  KDIE  value  emerges  for  the 
exothermic  portion  of  both  the  liquid  state  DNNC  and  RDX  thermochemical  decomposition 
processes. 

Table  6. 
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Comparative  KDIE  values  for  exothermic  liquid  state  DNNC  and  RDX  thermochemical 
decomposition  process 


Temp,  range 
(°C)  “ 

Reference 

Compound 

Analysis  method 

KDIE  value 

174-  194 

this  work 

DNNC/DNNC-d6 

isothennal  DSC 

1.65 

227  -  242 

12 

RDX/RDX-d6 

isothennal  DSC 

1.64 

212  only 

12 

RDX/RDX-d6 

isothennal 

QUVA" 

1.57 

206-216 

13 

RDX/RDX-d6 

isothermal  TGA 

1.56 

a  Quenched  Ultra-Violet  Absorbance  measurement 


For  comparison,  molecular  KDIE-based  IDSC  analyses  of  the  solid  state  DNNC 
thermochemical  decomposition  process  also  reveal  a  primary  1 .47  KDIE  and  a  rate¬ 
controlling  C-H  bond  rupture  for  the  exothennic  acceleratory  portion  of  the  decomposition 
process.  In  contrast  to  the  liquid  DNNC  decomposition  results,  however,  the  solid  state 
DNNC  induction  period  affords  a  secondary  1 .22  KDIE  value  suggesting  that  homolytic 
ring  C-N  ring  bond  cleavage  likely  constitutes  the  endothermic  rate-controlling  feature  prior 
to  transitioning  into  the  exothennic  acceleratory  phase  [2], 

The  homolytic  rate-controlling  chemical  bond  rupture,  and  resultant  atom  transfer  during 
the  decomposition  process  of  an  energetic  compound,  depends  upon  its  chemical  structure 
and  resultant  chemical  environment  [23].  As  seen  in  Fig.  9,  the  TNT  molecule  possesses  a 
chemical  structure  where  two  types  of  C-H  bonds  reside  in  vastly  different  chemical 
environments.  This  is  reflected  by  the  large  proton  NMR  shifts  between  the  C-H  bonds  of 
the  saturated  aliphatic  pendant  methyl  (-CH3)  group  and  the  aromatic  phenyl  C-H  bonds  (A 
=  6.32  ppm).  TNT  forms  mainly  condensed  phase  products  during  thermochemical 
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decomposition  [41,42]  via  a  rate-controlling  C-H  bond  rupture  at  the  pendant  methyl  group 
site  [10];  and  the  aromatic  C-H  bonds  remain  intact  when  fonning  these  condensed  phase 
decomposition  products  [9,10,41,42].  The  chemical  structure  of  the  DNNC,  where  all  C-H 
bonds  are  embodied  in  saturated  aliphatic  ring  methylene  groups,  represents  a  far  different 
structural  environment  from  TNT  (Fig.  9). 

Insert  Fig.  9 

The  rate-controlling  C-H  bond  rupture  found  in  the  exothermic  liquid  state  DNNC 
thermochemical  decomposition  process  only  can  emanate  from  the  saturated  methylene 
groups  in  the  DNNC  ring  structure  (Fig.  2);  but,  the  differences  in  their  chemical 
environment  are  much  more  subtle  than  with  TNT.  This  subtle  difference  in  DNNC 
molecular  structure,  and  resultant  chemical  environment,  is  illustrated  by  the  very  small 
NMR  proton  shift  between  the  two  C-2  methylene  protons  versus  the  four  C-4/C-6 
methylene  protons  (A=  0.83  ppm)  seen  in  Figure  9.  While  the  two  methylene  groups  at  the 
C-4  and  C-6  DNNC  ring  positions  are  chemically  and  structurally  equivalent,  the  methylene 
at  the  C-2  position  is  neither  chemically  nor  structurally  equivalent  with  the  C-4  and  C-6 
methylenes.  This  subtle  difference  in  the  DNNC  chemical  structure  raises  two  questions. 
Are  all  three  methylene  groups  in  the  DNNC  ring  structure  the  source  for  this  rate¬ 
controlling  C-H  bond  rupture  for  the  DNNC  thermochemical  decomposition  process,  like 
that  found  in  RDX  with  its  three  equivalent  methylenes?  Or,  secondly,  is  the  structural  and 
resultant  chemical  non-equivalence  of  the  DNNC  ring  structure  sufficiently  different  to 
restrict  this  C-H  bond  rupture  to  either  the  equivalent  C-4/C-6  methylene  sites,  or 
alternatively,  to  the  non-equivalent  nitrogen  sandwiched  C-2  site?  An  investigative 
mechanistic  IDSC  comparison  of  the  four  DNNC,  DNNC-d2,  DNNC-d4,  and  DNNC-d6 
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isotopomers  enters  into  the  DNNC  molecule  and  differentiates  between  these  chemically 
similar  DNNC  ring  methylenes.  In  doing  so,  it  determines  the  structural  site  from  which  C- 
H  bond  rupture  ultimately  regulates  the  DNNC  thennochemical  decomposition  process.  A 
similar  approach  was  used  for  the  nitroraromatic  HNS  (hexanitrostilbene)  compound  to 
analyze  the  more  chemically  diverse  aromatic  and  vinyl  C-H  bonds  contained  in  its 
molecular  structure,  but  no  primary  KDIE  or  resultant  rate-controlling  C-H  bond  rupture 
was  found  [43]. 

3.3  Structural  Kinetic  Deuterium  Isotope  Effect 

A  structural  KDIE  represents  an  approach  which  compares  different  internal  structural 
features  of  a  given  compound  by  using  both  the  unlabeled  compound  (DNNC)  and  the 
perdeuterio-labeled  compound  (DNNC-d6),  as  bracketed  baseline  data,  against  analogous, 
partially  labeled  isotopomers  (DNNC-d2,  DNNC-d/j).  There  is  a  9.6  kJ/mole  energy 
difference  between  the  C-H  bond  ground  state  vibrational  energy  and  that  of  the  analogous 
but  stronger  C-D  bond.  If  bond  energy  differences  between  the  non-equivalent  C-2 
methylene  and  the  equivalent  C-4/C-6  methylene  positions  in  the  unsymmetrical  DNNC 
molecule  are  greater  than  9.6  kJ/mole,  selectively  deuterium-labeled  DNNC-d2  and  DNNC- 
d4  isotopomers  could  reveal  which  methylene  position  provides  the  rate-controlling  C-H 
bond  rupture  for  the  DNNC  decomposition  process.  Consistent  with  the  KDIE  concept,  if 
the  DNNC-d2  and  DNNC-d4  isotopomers  seen  in  Fig.  10  were  synthesized  and  subjected  to 
IDSC  analyses,  the  methylene  position  providing  the  rate-controlling  C-H  bond  rupture 
could  be  identified.  Should  the  C-2  methylene  position  exclusively  provide  this  rate¬ 
controlling  C-H  bond  rupture,  then  the  partially  labeled  DNNC-d2  isotopic  isomer  would 
give  the  same  slower  rate  constant  (kd2)  as  DNNC-d6  (kd6),  and  the  partially  deuterium 
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labeled  DNNC-d4  would  give  the  same  rate  constant  (kd4)  as  DNNC  (k).  Conversely,  if  the 
chemically  equivalent  C-4/C-6  methylenes  are  the  source  of  the  C-H  rate-controlling  bond 
cleavage,  DNNC-d2  would  have  the  same  rate  constant  (kd2)  as  DNNC  (k),  and  DNNC-d4 
(kd4)  should  give  the  same  rate  constant  as  does  the  slower  one  exhibited  by  DNNC-d6  (kd(-,). 

Insert  Fig.  10 

The  partially  labeled  isotopic  isomers,  DNNC-d2  and  DNNC-d4,  were  synthesized 
(Experimental  Section)  where  discreet  perdeuterio  -CD2-  methylene  groups  are  introduced 
into  the  DNNC  molecular  structure.  To  the  nearest  whole  number,  proton  NMR  integration 
verified  that  each  contained  the  intended  two  and  four  deuterium  atoms,  respectively. 
However,  for  both  the  synthesized  DNNC-d?  and  DNNC-d4  compounds,  proton  NMR 
integration  for  the  specific  C-2  and  C4/C-6  methylenes  showed  a  1:2  integration  ratio 
indicating  that  random  -CD2-  methylene  scrambling  into  all  three  DNNC  ring  positions 
resulted  with  both  isotopomers.  Under  the  alkaline  reaction  conditions,  it  was  found  that  all 
three  C-2,  C-4,  and  C-6  positions  in  the  DNNC  ring  structure  are  subject  to  a 
formylation/defonnylation  equilibrium  involving  intermixing  fonnaldehyde  molecules 
which  are  generated  from  the  ADIOL  reactant  and  available  from  the  additional 
formaldehyde  reactant.  Unfortunately,  this  equilibrium  exchange  occurred  faster  than  did 
ring  self-assembly  (Fig.  1 1). 

Insert  Fig.  1 1 

Mass  spectral  analysis  of  the  synthesized  DNNC-d2  and  DNNC-d4  samples  revealed  that 
each  contained  a  mixture  of  all  four  isotopomers,  DNNC,  DNNC-d2,  DNNC-d4,  and  DNNC- 
d6,  in  percentages  that  mirrored  the  intended  pure,  selectively-positioned  isotopic  isomers 
DNNC-d2  and  DNNC-d4  seen  in  Fig.  10.  The  DNNC-d2  isotopomer  predominated  in  the 
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mixed,  synthesized  DNNC-d?  sample,  and  the  DNNC-d4  isotopomer  was  predominant  in  the 
mixed,  synthesized  DNNC-d4  sample.  The  net  result  of  this  random  -CD2-  methylene  ring 
incorporation,  and  the  resultant  percent  mixtures  of  all  four  isotopic  isomers  in  the 
synthesized  DNNC-d?  and  DNNC-d4  samples,  produced  bulk  materials  where  the  DNNC-d? 
behaved  as  if  it  were  one-third  deuterium  atom  labeled  at  each  methylene  (C-2,  C-4,  &  C-6) 
ring  position,  and  the  DNNC-d4  bulk  material  behaved  as  if  it  were  two-thirds  deuterium 
labeled  at  each  methylene  position  (Fig.  12).  This  serendipitous  result  necessitated  using  a 
modified  and  slightly  more  complex  approach  to  detennine  whether  or  not  a  specific 
methylene  group  in  the  DNNC  chemical  structure  is  responsible  for  the  rate-controlling  C-H 
bond  rupture 

Initially,  it  was  surmised  that  such  random  isotopic  methylene  group  scrambling 
possibly  explained  the  partial  rate  reductions  seen  as  one  progresses  from  DNNC  to  DNNC- 
d2,  DNNC-d4,  and  finally  to  DNNC-d6  (Table  7)  without  revealing  any  mechanistic 
differentiation  between  the  two  dissimilar  DNNC  methylene  groups  [44].  An  evolving 
mechanistic  analysis,  however,  suggested  otherwise. 

Random  isotopic  methylene  scrambling  could  only  give  this  result  if  the  DNNC  molecule 
predominantly  decomposes  with  the  ring  structure  immediately  unraveling  to  form  one 
neutral  CH2=CN02)2,  and  two  CH2=NN02  fragments,  or  even  lower  molecular  weight 
species,  prior  to  any  rate-controlling  C-H  bond  rupture  occurring.  If  this  were  the  case,  a 
molecular  KDIE  reflecting  a  secondary  value  should  have  been  observed.  Observation  of  a 

Table  T. 

IDSC  rate  constant  data  reduction  for  the  DNNC  and  its  deuterated  isotopomers 
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Kinetic  data  DNNC 


DNNC-d2 


DNNC-d4 


DNNC-d6 


Rate  constant  7.79  x  10'3  6.84  x  10'3  5.92  x  10'3  4.68  x  10  3 

(kins')  +0.25  ±0.30  ±0.57  ±0.61 

Interval  (s'1)  0.95  x  10'3  0.92  x  10'3  1.24  x  10'3 

Percent  k  0%  31%  60%  100% 

reduction 


Accel,  phase  93  89  88  76 

linearity  (%) 

No.  runs  6  5  5  5 

a  All  data  variances  computed  at  the  99%  Student  t  level 


primary  KDIE  value  strongly  suggests  this  is  not  the  case  and  that  the  skeletal  ring  structure 
primarily  stays  intact  through  the  initial  and  subsequent  rate-controlling  C-H  bond  rupture. 

Secondly,  past  condensed  phase  product  isolation  studies  with  TNT  [41,42]  and  TATB 
[19]  show  intact  ring  structures  where  decomposition  reactions  occur  only  with  the  pendant 
methyl,  amino,  and  nitro  groups  and  include  those  consistent  with  the  KDIE-detennined  rate 
controlling  C-H  (TNT)  [10]  and  N-H  (TATB)  [3]  bond  ruptures.  Ring  structures  are  left 
intact,  and  initial  reactions  with  nitramines  appear  to  involve  the  pendant  NCE  group  [5,8] 
and  subsequent  rate-controlling  ring  hydrogen  atom  rupture  by  NO?  radical  abstraction  [8]. 

Thirdly,  sub-initiation  impact  and  shock  studies,  where  the  charred  inner  surfaces  of  hot 
spots  caused  by  a  rapid  thermochemical  decomposition  process,  also  reveal  condensed  phase 
intennediate  condensed  phase  products  with  intact  ring  structures  for  TNT  [22],  TATB  [21], 
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and  RDX  [20,44].  In  the  case  of  structurally  related  RDX,  sub-initiation  by  impact  displays 
trigger  linkage  or  first  reaction  N-0  bond  rupture  to  form  the  analogous  trinitroso 
intennediate  compound  [45].  More  drastic  sub-initiated  shock  reveals  a  melamine  type  of 
product  where  C-H  bond  rupture  [20],  the  rate-controlling  step  reflected  by  KDIE-based 
thermochemical  decomposition  investigations  [1 1,12],  apparently  follows  the  initial  N-0 
bond  rupture,  and  again,  leaves  the  skeletal  ring  structure  intact  as  discussed  in  the 
Introduction  section. 

A  primary  molecular  KDIE,  involving  a  rate-controlling  C-H  bond  rupture,  occurs  early 
during  the  thermochemical  decomposition  process  and  involves  an  intennediate  DNNC 
radical  with  an  intact  skeletal  ring  structure.  Literature  precedence  also  suggests  that  the 
skeletal  ring  structure  of  cyclic  energetic  materials  is  attacked  late  during  the 
thennochemical  decomposition  process  after  the  rate-controlling  bond  rupture  occurs. 
Therefore,  the  consistent  decrease  in  DNNC  decomposition  rate  with  increased  deuterium 
substitution  in  Table  7  cannot  be  a  random  occurrence  under  these  documented  mechanistic 
circumstances.  A  more  complex  modified  evaluative  approach  is  required  for  mechanistic 
consistency  [30], 

Insert  Fig.  12 

This  modified  approach  invokes  two  reasonable  assumptions.  First,  when  possible,  the 
homolytic,  covalent  C-H  bond  rupture,  and  H  atom  transfer,  will  always  occur  in  preference 
to  that  of  the  stronger  C-D  covalent  bond.  Secondly,  should  the  rate-controlling  step 
involve  and  be  forced  to  proceed  at  one  specific  DNNC  methylene  site,  and  this  methylene 
site  is  deuterium  labeled,  then  C-D  bond  rupture  will  occur  with  a  concomitant  decrease  in 
the  decomposition  rate  and  associated  rate  constant. 
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Fig.  12  represents  a  valid  model  with  which  to  evaluate  kinetics  data  obtained  with 
DNNC  and  its  three  deuterium-labeled  isotopomers,  DNNC-d2,  DNNC-d/t,  and  DNNC-d6. 
Fig.  12  (top)  displays  three  -CD2-  randomly  labeled  DNNC-d2  molecules,  which  combined, 
mirror  the  synthesized  DNNC-d2  isotopic  isomer  that  behaves  as  if  it  were  one-third 
deuterium  labeled  at  each  methylene  position.  If  the  rate-controlling  C-H  bond  rupture 
emanates  from  the  C-2  methylene  site,  slower  C-D  bond  cleavage  would  occur  only  one  out 
of  three  times,  and  C-H  bond  rupture  would  occur  twice  as  often.  Thus,  DNNC-d2  should 
give  a  one  third  (33%)  rate  constant  reduction  between  the  faster  unlabeled  DNNC  rate 
constant  and  the  relative  100%  rate  constant  reduction  of  DNNC-d6.  Alternatively,  if  the 
equivalent  C-4/C-6  methylene  positions  are  the  source  for  the  rate-controlling  C-H  bond 
rupture,  a  slower  C-D  bond  rupture,  and  concomitant  rate  constant  reduction  would  not 
occur  with  the  DNNC-d2  isotopic  isomer  because  a  preferred  C-H  bond  rupture  is  always 
available.  Fig.  12  (bottom)  displays  three  -CD2-  randomly  labeled  DNNC-d4  molecules, 
which  combined,  mirror  the  synthesized  DNNC-d4  isotopic  isomer  that  behaves  as  if  it  were 
two-thirds  deuterium  labeled  at  each  methylene  position.  With  the  synthesized  DNNC-cU 
isotopic  isomer,  a  C-H  rate-controlling  bond  rupture  relegated  to  the  C-2  methylene  would 
effect  a  two-thirds  (67%)  rate  constant  reduction  since  this  methylene  site  is  deuterated  two 
out  of  three  times.  If  the  rate-controlling  step  were  provided  by  the  C-4/C-6  sites,  DNNC-d4 
would  give  a  one  third  (33%)  rate  constant  reduction.  This  is  summarized  in  Table  8. 


Table  8. 

Expected  rate  constant  reduction  of  methylene  sites  relative  to  DNNC/DNNC-d6  rate 
constant  difference 
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Isotopomer  compound 


RCS  at  C-2 

Rate  constant  reduction 


RCS  at  C-4/C-6 

Rate  constant  reduction 


DNNC-d2 

33% 

0% 

DNNC-d4 

67% 

33% 

DNNC-d6 

100% 

100% 

RCS  =  rate-controlling  step  (C-H  bond  rupture) 

The  DNNC,  DNNC-dj,  DNNC-d4,  and  DNNC-d6  thermochemical  decomposition 
processes  were  analyzed  by  IDSC  at  186  °C,  and  a  systematic  rate  reduction  was  observed. 
Fig.  13  shows  an  individual  IDSC  run  for  each  isotopic  isomer  where  the  decomposition  rate 
progressively  slows  in  the  order:  (fastest)  DNNC  >  DNNC-d2  >  DNNC-d4  >  DNNC-d6 
(slowest).  The  fact  that,  both  the  DNNC-d2  and  DNNC-d4  rate  constants  are  reduced 
significantly  (Table  8)  during  the  exothennic  acceleratory  phase,  immediately  identifies  the 
C-2  methylene  position  as  being  the  structural  site  for  the  rate-controlling  C-H  bond  rupture 
for  this  portion  of  the  thermochemical  decomposition  process. 

Insert  Fig.  13 

Careful  rate  constant  data  analysis  of  the  experimental  structural  KDIE  data  in  Table  8 
shows  that,  within  the  total  rate  difference  between  DNNC  and  fully  deuterium-labeled 
DNNC-d6,  DNNC-d2  decomposes  about  one-third  (31%)  more  slowly  than  DNNC,  and 
DNNC-d4  nearly  two-thirds  (60%)  more  slowly  (Table  8).  That  is,  the  structural  KDIE  rate 
reduction  percent  for  DNNC-d2  equals  [kDNNC  -  kDNNC-d2  divided  by  kDNNC  -  kDNNC-d6]  times 
100,  and  the  structural  KDIE  rate  reduction  percent  for  DNNC-d4  equals  [kDNNC  -  kDNNC-d4 
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divided  by  kDNNC  -  kDNNC-d6]  times  100.  This  data  very  closely  matches  the  calculated 
criteria  expected  if  the  non-equivalent  C-2  methylene  position  in  the  DNNC  molecule  is  the 
origin  for  the  rate-controlling  C-H  bond  rupture  site,  and  further  experimental  data  analysis 
strengthens  this  conclusion. 

Rate  constant  interval  differences  in  Table  8  reveal  nearly  identical  interval  reductions 
between  the  DNNC  and  DNNC-d2  rate  constant  values,  and  between  the  DNNC-d2  and 
DNNC-d4rate  constants;  however,  a  much  larger  1.24  X  10'  s'  interval  reduction  appears 
between  the  DNNC-d4  and  DNNC-d6  rate  constants.  Fig.  14  displays  this  data  in  graph  plot 
form  and  further  suggests  the  experimental  DNNC-d6  rate  constant,  while  quite  reasonable, 
might  actually  be  slightly  lower  than  what  is  consistent  with  the  other  DNNC,  DNNC-d2, 
and  DNNC-d4  isotopomers.  Plotting  the  number  of  hydrogen  atoms  (x-axis)  in  each  DNNC 
isotopomer  [46],  versus  the  measured  rate  constant  (y-axis)  in  Figure  14,  affords  an 
excellent  linear  least  squares  fit  (R  =  0.994),  but  Figure  14  also  gives  a  linear-intercept 
value  for  DNNC-d6  (  0  H- Atoms)  [46]  which  suggests  a  rate  constant  equal  to  4.77  x  10'  s' 

1  3  1 

,  rather  than  the  lower  4.68  x  10'  s'  ,  would  be  more  reasonable  [47]. 

Insert  Fig.  14 

By  graphing  only  the  first  three  points  corresponding  to  DNNC,  DNNC-d2,  and  DNNC- 
d4  in  an  analogous  plot  (Fig.  15),  a  near  perfect  linear  least  squares  fit  (R  =  0.9999  or 
1 .000)  results  and  gives  a  y-intercept  value  at  the  DNNC-d6  (0  H-atoms)  isotopomer  that 
defines  a  4.98  x  10'  s'  revised 

Insert  Fig.  15 
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3  1 

rate  constant  value.  When  this  revised  rate  constant  (4.98  x  10'  s'  )  is  substituted  for  the 
4.68  x  10"  s'  experimental  data  rate  constant  in  Table  9,  all  three  rate  constant  reduction 
interval  differences  match  much  more  closely  with  the  re-computed  DNNC-d4  to  DNNC-d6 
interval  falling  between  the  two  other  intervals.  Most  importantly,  a  34%  rate  constant 
reduction,  versus  calculated  33%,  results  for  DNNC-d2,  and  the  expected  67%  rate  reduction 
surfaces  for  DNNC-d4.  Thus,  the  rate-controlling  C-H  bond  rupture,  which  regulates  the 
energy  release  rate  of  the  exothermic  acceleratory  phase,  originates  at  the  non-equivalent  C- 
2  methylene  site  sandwiched  between  the  two  DNNC  ring  nitrogen  atoms. 

Table  9. 

Revised  DNNC-d6  rate  constant  data  for  DNNC  and  its  deuterated  isotopomers 


Kinetic  data 

DNNC 

DNNC-d2 

DNNC-d4 

DNNC-d6 

Rate  constant 

7.79  x  10'3 

6.84  x  10'3 

5.92  x  10'3 

4.98  x  10  3 

(k  in  s'1) 

Interval  (s'1) 

0.95  x  10'3 

0.92  x  10'3 

0.94  x  10'3 

Percent  k 
reduction 

0% 

34% 

67% 

100% 

Next,  the  induction  period  data  is  similarly  evaluated  to  determine  if  this  mechanistic 
feature  also  regulates  this  earlier  portion  of  the  DNNC  thermochemical  decomposition 
process.  Induction  period  times  also  represent  rate  measurement  data  and  have  an  inversely 
proportional  relationship  with  corresponding  reaction  rate  constants  which  cannot  be 
directly  measured  by  IDSC  in  this  flat  portion  of  the  thermogram  [10].  A  fast  endothermic 
reaction  rate  (large  rate  constant),  occurring  during  the  induction  period,  provides  a  short 
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(small)  induction  time.  Therefore,  induction  times  should  increase  as  the  amount  of 
deuterium  labeled  methylene  groups  increase  in  the  DNNC  isotopomer  by  the  proportionate 
percentages  shown  in  parentheses:  DNNC  (0%  baseline  time)  <  DNNC-d2  (33%)  <  DNNC- 
d4  (67%)  <  DNNC-d6  (100%  maximum  time). 

Table  10  lists  the  averaged  induction  time  (h)  obtained  at  186  °C  for  DNNC  and  each  of 
its  three  deuterium  labeled  isotopomers.  The  increased  induction  times  for  both  DNNC-d? 
and  DNNC-d4,  compared  to  DNNC  (Fig.  16),  again  suggests  that  the  rate-controlling  C-H 
bond  rupture  emanates  from  the  C-2  methylene  site  of  the  DNNC  ring  structure.  Should  this 
rate-controlling  mechanistic  feature  be  generated  from  the  equivalent  C-4  and  C-6 
methylenes,  no  induction  time  decrease  for  the  DNNC-d2  isotopomer,  compared  to  DNNC, 
should  result  (Fig.  12).  Induction  time  interval  differences  between  DNNC/DNNC-d2, 
DNNC-d2/DNNC-d4,  and  DNNC-d4/DNNC-d6  do  not  follow  the  expected  regular 
magnitude  increase  which  should  give  a  33%  induction  time  increase  for  DNNC-d2  and  a 
67%  increase  for  DNNC-d4.  Rather,  induction  time  increases  are  29  %  for  DNNC-d2  and 
46%  for  DNNC-d4  using  the  initial  experimental  interval  data  obtained  (Table  10). 


Table  10.fl 

Induction  time  data  for  DNNC  and  its  partially  deuterated  isotopomers  at  186  °C 
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Induction 
period  data 


DNNC  Interval  DNNC-d2  Interval  DNNC-d4  Interval  DNNC-d6 


No.  IDSC  6  5  5  5 

runs 


Ave.  91.8  ±23.0  115.7±13.8  129.2  ±  25.1  173.9  ±34.1 

induction 

time  (s) 

Interval  (s)  23.9  13.5  44.7 

Percent  0%  29%  46%  100% 

reduction 

a  All  data  variances  computed  at  the  99  percent  Student  t  level 


Again,  like  the  exothermic  acceleratory  phase  data,  an  extended  treatment  of  the 
induction  time  data  clearly  reflects  the  expected  trend.  A  graphic  trend  line  plot  of 
experimental  induction  times  versus  the  number  of  hydrogen  atoms  in  DNNC,  DNNC-d2, 
DNNC-d.4,  and  DNNC-d6  (Table  10),  only  gives  an  R2  =  0.946  linear  least  squares  fit  (Fig. 
16);  but,  comparison  of  a  polynomial  curved  trend  line  and  close  examination  of  the  linear 
trend  line  and  interval  differences  from  Table  10,  suggest  the  induction  time  for  DNNC-d4 
may  be  low  compared  to  the  other  three  isotopomers. 

Insert  Fig.  16 

In  contrast  to  Fig.  16,  Fig.  17  plots  a  least  squares  trend  line  fit  for  the  induction  times  of 
DNNC,  DNNC-d2,  and  DNNC-d6  only,  where  the  DNNC-d4  129.2  s  induction  time  data  is 
excluded,  and  affords  an  excellent  R~  =  0.998  linear  fit  with  a  13.81  slope  value.  A  similar 
plot  that  includes  the  DNNC,  DNNC-d2,  and  DNNC-d4  induction  times,  but  excludes  the 
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2 

DNNC-d6  173.9  s  induction  time  (Fig.  18),  gives  a  much  lower  R  =  0.976  fit,  and  verifies 
that  the  abnormally  low  induction  time  comes  from  DNNC-d4  rather  than  DNNC-d6. 

Insert  Fig.  17  &  18 

Averaging  the  sum  of  the  three  induction  time  interval  values  found  in  Table  10,  an 
interval  of  27.4  s  results  over  the  entire  DNNC  to  DNNC-d6  isotopomer  range.  Substituting 
this  averaged  27.4  s  interval  difference  for  the  questionable  13.5  second  interval  difference 
between  DNNC-d?  and  DNNC-d4,  and  adding  the  27.4  s  to  the  1 15.7  s  induction  time  of 
DNNC-di,  a  143.1  s  induction  time  for  DNNC-d4 results  which  nearly  matches  the  upper 
99%  Student  t  variance  for  this  questionable  129.2  s  induction  period  value.  The  calculated 
27.4  s  induction  time  interval,  and  resultant  calculated  143.1  s  induction  time  for  DNNC-d4, 
along  with  a  revised  30.8  s  interval  difference  in  place  of  44.7  s  between  DNNC-d4  and 
DNNC-d6  to  maintain  a  173.9  s  induction  time  for  DNNC-d6,  are  substituted  for  the 
corresponding  data  in  Table  10  to  give  Table  11.  Using  the  revised  data  from  Table  11,  the 
DNNC-d2  induction  time  increase  remains  at  29%  for  DNNC-d2,  but  becomes  a  much  more 
reasonable  62%  for  DNNC-d4.  The  29%  DNNC-d2  induction  time  increase  is  reasonably 
close  to  the  expected  33%  value,  and  the  corrected  62%  induction  time  increase  for  DNNC- 
d4  (Table  1 1),  versus  46%  from  Table  10,  compares  more  favorably  with  the  calculated  67% 
induction  time  increase  expected. 


Table  11. 

Revised  induction  time  data  for  DNNC  and  its  partially  deuterated  isotopic  isomers 


Induction  DNNC  Interva  DNNC-d2  Interva  DNNC-d4  Interva  DNNC-d6 
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period  data 


1 


1 


1 


No.  IDSC  6  5  5  5 

Runs 

Ave.  91.8  115.7  143.1  173.9 

induction 

time  (sec) 

Interval  23.9  27.4  30.8 

(sec) 

Percent  0%  29%  62%  100% 

reduction 

A  graphic  plot  of  the  revised  induction  time  data  in  Table  1 1  for  DNNC,  DNNC-d2, 
DNNC-d.4,  and  DNNC-d6  generates  a  linear  trend  line  with  a  nearly  identical  linear  R“  = 
0.997  fit  and  13.68  slope  (Fig.  19)  as  that  seen  with  the  linear  R2  =  0.998  value  and  13.81 
slope  found  in  Fig.  16  where  the  suspect  DNNC-d4  induction  time  is  excluded.  Allowing 
that  visually  interpreted  induction  time  data  measurements  are  not  as  exact  as  acceleratory 
phase  rate  constants,  this  induction  time  data  treatment  further  confirms  that  C-H  bond 
rupture  at  the  DNNC  structural  C-2  methylene  position  also  controls  the  endothermic 
induction  period  portion  of  the  liquid  state  DNNC  thermochemical  decomposition  process. 

Insert  Fig.  19 


4.  Conclusions 
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Liquid  state  thermochemical  decomposition  of  DNNC,  conducted  between  174  to  194 
°C,  is  an  autocatalytic  process.  During  this  process  homolytic  covalent  C-H  bond  rupture, 
and  resultant  hydrogen  atom  transfer,  constitute  a  key  mechanistic  feature  which  detennines 
the  endothermic  initiation  and  subsequent  exothermic  propagating  energy  release  rates. 

Molecular  KDIE  data,  taken  with  un-labeled  DNNC  and  perdeuterio-labeled  DNNC-d6, 
displays  a  primary  KDIE  value  for  both  the  endothermic  induction  period  and  subsequent 
autocatalytic  exothennic  acceleratory  phase  (Fig.  20,  top).  The  initial  endothermic 
induction  period  likely  generates  a  catalytic  species  via  a  rate-controlling  homolytic  C-H 
bond  rupture  (KDIE  =  1.76),  which  upon  reaching  a  threshold  concentration,  initiates  an 
exothermic  autocatalytic  decomposition  process.  The  exothermic  decomposition  process 
then  follows  a  rate  that  also  is  controlled  by  a  homolytic  C-H  bond  rupture  (KDIE  =  1 .65) 
coming  from  one  of  the  methylene  groups  comprising  the  DNNC  ring  structure.  A 
significantly  higher  energy  of  activation  is  required  to  initiate  the  DNNC  endothermic 
induction  period  (170.0  kJ/mol)  than  what  is  needed  for  the  corresponding  exothermic 
acceleratory  portion  of  the  respective  decomposition  process  (104.2  kJ/mol).  The  65.8 
kJ/mol  energy  difference,  coupled  with  its  autocatalytic  rate  behavior  supports  this  catalytic 
mechanistic  interpretation  for  the  liquid  state  DNNC  thermochemical  decomposition 
process.  This  same  mechanistic  trend  previously  was  reported  for  the  solid  state  DNNC 
thennochemical  decomposition  process  (Table  3)  [2], 

Insert  Fig.  20 

Structural  KDIE  data,  obtained  from  comparing  DNNC  and  perdeuterio-labeled 
DNNC-d6  with  partially  deuterium  labeled  isotopomers  DNNC-d?  and  DNNC-dzi,  enters  into 
the  DNNC  molecular  ring  structure  and  reveals  that  the  rate  controlling  C-H  bond  rupture 
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originates  at  the  methylene  (-CH2-)  site  at  the  C-2  position  of  the  DNNC  molecule  which  is 
sandwiched  between  the  two  ring  nitrogen  atoms,  or  originates  from  a  structural  fragment 
that  is  generated  from  this  structural  site  (Fig.  20,  bottom).  The  chemically  and  structurally 
equivalent  methylene  groups  at  the  C-4  and  C-6  positions  of  the  DNNC  molecule  reside  in  a 
chemical  environment  which  is  subtly  different  from  the  third  C-2  methylene  group.  The 
subtle  energy  difference  of  the  chemical  environments  between  the  C-2  methylene  site  and 
the  C-4/C-6  equivalent  methylene  sites,  or  the  corresponding  structural  fragments  that 
provide  the  rate-controlling  C-H  bond  rupture,  must  exceed  the  9.6  kJ/mol  ground  state 
vibrational  energy  difference  between  the  C-H  and  C-D  bond.  If  this  were  not  the  case,  the 
rate-controlling  C-H  bond  rupture,  and  hydrogen  atom  transfer,  would  occur  randomly  at  all 
three  DNNC  methylene  positions.  The  C-4/C-6  methylenes  contribute  little  or  no 
significant  role  in  determining  the  overall  liquid  state  DNNC  thermochemical 
decomposition  and  energy  release  rates. 

Apparently,  the  liquid  state  DNNC  induction  period  requires  a  higher  activation  energy 
(170  kJ/mol)  than  does  the  solid  state  (128  kJ/mol)  [2]  to  initiate  the  subsequent  exothermic 
portion  of  the  thermochemical  decomposition  process.  The  liquid  state  (104  kJ/mol)  and 
solid  state  (105  kJ/mol)  [2]  DNNC  exothermic  acceleratory  phases  have  a  similar  activation 
energy  once  energy  release  commences.  While  homolytic  C-H  bond  rupture  determines 
both  the  endothermic  induction  time  duration  and  the  exothermic  acceleratory  phase  energy 
release  rate  in  liquid  state  DNNC  thermochemical  decomposition,  rate-controlling  homolytic 
C-H  bond  dissociation  only  is  found  during  the  acceleratory  phase  of  the  solid  state 
thermochemical  decomposition  process.  Solid  state  DNNC  endothennic  induction  period 
kinetics  appear  to  be  regulated  by  a  different  rate-controlling  bond  rupture  as  previously 
reported  [2]. 
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Figure  1.  General  thermochemical  decomposition  process  schematic  representation  and  the 
KDIE  measurement  context  [13]. 
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Figure  2.  Unique  hybrid  chemical  structure  of  the  DNNC  molecule. 


Distribution  A.  Approved  for  public  release;  distribution  unlimited. 


43 


Figure  3.  Individual  IDSC  thennograms  of  DNNC  from  174  °C  to  194  °C  taken  at  4  °C 
intervals. 
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da/dt 


Figure  4.  (top)  DNNC  isothermal  DSC  thermogram  (Run  D6H26)  illustrating  the  portion 
of  the  exothennic  acceleratory  phase  that  fits  autocatalytic  kinetics  following  the 
endothermic  induction  period  (t,).  (bottom)  Corresponding  autocatalytic  rate  plot  for 
DNNC  (Run  D6H26)  illustrating  the  linear  kinetic  fit  of  the  exothermic  acceleratory  phase 
and  slope-determined  rate  constant  (k). 
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Figure  5.  Kinetic  data  plot  providing  the  activation  energy  for  the  DNNC  endothermic 
induction  period. 
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Figure  6.  Kinetic  data  plot  for  activation  energy  of  the  DNNC  &  DNNC-d6  exothermic 
acceleratory  phase. 
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Figure  7.  Individual  DNNC  and  DNNC-d6  IDSC  thermograms  compared  at  174  °C. 
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Figure  8.  Possible  KDIE  types  from  C-H  (1°  KDIE)  bond  rupture  and  C-N,  N-N02,  C-NO? 
(2°  KDIE)  potentially  present  with  the  DNNC  molecular  structure. 


Distribution  A.  Approved  for  public  release;  distribution  unlimited. 


49 


I 


Figure  9.  Comparative  proton  NMR  chemical  shift  magnitudes  of  TNT  C-H  bonding  (A= 
6.32  ppm)  versus  DNNC  ring  methylene  C-H  bonding  (A=  0.83  ppm). 
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Figure  10.  Theoretical  KDIE  results  with  specific  partially  deuterium-labeled  DNNC-di 
and  DNNC-d4  isotopomers. 
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Figure  11.  General  synthesis  scheme  used  for  ADIOL,  ADIOL-d4  and  the  four  DNNC 
isotopomers. 
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Figure  12.  Overall  kinetic  behavior  of  random  partially  deuterium-labeled  DNNC-d2  and 
DNNC-d4  isotopomer  bulk  samples. 
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HEAT  EVOLVED  (mW) 


Figure  13.  Comparison  of  individual  DNNC,  DNNC-d2,  DNNC-d4,  and  DNNC-d6 
isotopomer  thermograms  at  186  °C. 
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RATE  (SEC-1) 


RATE  CHANGE  VS.  DNNC  D-ATOM 
INCORPORATION 
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Figure  14.  Experimental  DNNC  isotopomer  rate  constant  versus  number  of  methylene 
hydrogen  atoms  in  DNNC  molecular  structure. 
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RATE  (SEC-1) 
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INCORPORATION 
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Figure  15.  Experimental  DNNC  isotopomer  rate  constant  versus  number  of  methylene 
hydrogen  atoms  in  DNNC  molecular  structure  with  DNNC-d6  (H  =  0  atoms)  isotopomer  rate 
constant  not  plotted. 
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Figure  16.  Experimental  DNNC  isotopomer  induction  period  times  versus  number  of 
methylene  deuterium  atoms  in  DNNC  molecular  structure. 


Distribution  A.  Approved  for  public  release;  distribution  unlimited. 


57 


TIME  CHANGE  VS.  D-ATOM  INCORPORATION 
(RAW  DATA:  EXCLUDING  DNNC-d4) 
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Figure  17.  Experimental  DNNC  isotopomer  induction  period  time  versus  number  of 
methylene  deuterium  atoms  in  DNNC  molecular  structure  with  DNNC-d4  isotopomer 
induction  time  not  plotted. 
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TIME  CHANGE  VS.  D-ATOM  INCORPORATION 
(RAW  DATA:  EXCLUDING  DNNC-d6) 
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Figure  18.  Experimental  DNNC  isotopomer  induction  period  time  versus  number  of 
methylene  deuterium  atoms  in  DNNC  molecular  structure  with  DNNC-d6  isotopomer 
induction  time  not  plotted. 


Distribution  A.  Approved  for  public  release;  distribution  unlimited. 


59 


Figure  19.  Experimental  DNNC  isotopomer  induction  period  times  versus  number  of 
methylene  deuterium  atoms  in  DNNC  molecular  structure  with  DNNC-d4  isotopomer 
modified  induction  time  plotted. 
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Figure  20.  Graphic  conclusion  for  liquid  state  DNNC  thermochemical  decomposition 
findings. 
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